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Abstract 

We compute the dominant two- loop electroweak correction, of 0{GfM^), to 
the partial width of the decay of an intermediate-mass Higgs boson into a pair of 
photons. We use the asymptotic-expansion technique in order to extract the leading 
dependence on the top-quark mass plus four expansion terms that describe the 
dependence on the W- and Higgs-boson masses. This correction reduces the Born 
result by approximately 2.5%. As a by-product of our analysis, we also recover the 
0{GFMf ) correction to the partial width of the Higgs-boson decay to two gluon 
jets. 
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1 Introduction 



Among the main tasks of the current experiments at the Fermilab Tevatron and the 
future experiments at the CERN Large Hadron Colhder (LHC) is the search for the Higgs 
boson, which is the only missing particle in the standard model (SM). The electroweak 
precision data mainly collected at CERN LEP and SLAC SLC in combination with the 
direct top-quark mass measurement at the Tevatron favour a light Higgs boson with mass 
Mh = 113142 GeV with an upper bound of about 237 GeV at the 95% confidence level 
The allowed mass range is compatible with the so-called intermediate mass range, defined 
by Mw < Mh < 2M\y. In this mass range, the decay into two photons represents one of 
the most useful detection modes at hadron colliders. 

Since there is no direct coupling of the Higgs boson to photons, the process H — > 
77 is loop-induced. In the limit of vanishing bottom-quark mass, one distinguishes at 
lowest order the contributions from virtual top quarks and W bosons, where, in covariant 
gauge, the latter are accompanied by charged Goldstone bosons (0) and Faddeev-Popov 
ghosts {u). Some sample Feynman diagrams are depicted in Fig. ^ The corresponding 
contributions have been evaluated for the first time in Ref. j2] (for reviews, see also 
Ref. 13). QCD corrections, which only affect the top-quark diagrams are known at the 
two- |lj and three-loop j3] orders. Recently, also the two-loop electroweak correction 
induced by light-fermion loops has been evaluated pj. In this paper, we compute the two- 
loop electroweak correction that is enhanced by GpM^. For this purpose, we consider 
the (formal) hierarchy ^ {2M]y)^ ^ Mjj and apply the method of asymptotic 
expansion [7], which allows us to also evaluate four expansion terms in the ratio tw = 
Mh/{2MwY aside from the leading term in M"} . 

Due to electromagnetic gauge invariance, the amputated transition-matrix element of 
if — s> 77 possesses the structure 

T^' = {qi-q2g'''-qU'2)A (1) 

where /i and v are the Lorentz indices of the external photons with four-momenta qi and 
g2, respectively. Thus, the decay rate of the Higgs boson into two photons is given by 

r(i/-77) = ^I^P. (2) 
The form factor A is evaluated in perturbation theory as 

A = AfHA^S)+^m + ---. (3) 

where Af^ and Ay) denote the one-loop contributions induced by virtual top quarks 
and W bosons, respectively, A^^ stands for the two-loop electroweak correction involving 
virtual top quarks, and the ellipsis represents the residual one- and two-loop contributions 
as well as all contributions involving more than two loops. 

In the practical calculation, it is convenient to project out the scalar amplitudes that 
multiply the basic Lorentz tensors g^" , qiQ^^ gi5'2- The corresponding projectors 
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Figure 1: Sample Feynman diagrams contributing in leading order to the process H — >■ 77. 



can be obtained by in turn contracting T^'^ with these Lorentz tensors and solving the 
resulting system of linear equations. We separately project out the coefficients of the 
tensors qrq2g^'^ and qiq2i and thus have a strong check on our calculation. Furthermore, 
we adopt a general gauge in our calculation and verify that the gauge parameter drops 
out in the final result. For simplicity, the element Vtb of the Cabibbo-Kobayashi-Maskawa 
quark mixing matrix is set to unity, so that the quarks of the third fermion generation 
decouple from those of the first two, which they actually do to very good approximation 

n 

Our paper is organized as follows. In Section 12 we illustrate the usefulness of the 
asymptotic-expansion technique by redoing the one-loop calculation. In Section El we 
discuss the two-loop calculation. SectionlUcontains the discussion of the numerical results. 
We conclude with a summary in Section El 

2 One-loop results and counterterm contribution 

The application of the asymptotic-expansion technique to the one-loop diagrams, some of 
which are shown in Fig. ^ leads to a naive Taylor expansion in the external momenta qi 
and q2- Nevertheless, we use already here a completely automated set-up, which consists in 
the successive use of the computer programs QGRAF 9j, q2e [TUl, exp and MAT AD [T2]. 
First, QGRAF is used to generate the Feynman diagrams. Its output is then rewritten by 
q2e to be understandable by exp. In the two-loop case, the latter performs the asymptotic 
expansion and generates the relevant subdiagrams according to the rules of the so-called 
hard- mass procedure j7j. Form [13j files are generated, which can be read by MAT AD [T^ . 
which performs the very calculation of the diagrams. 

The analytic expression for the Born result can be found in Refs. j21 For com- 
pleteness, we list it here in closed form and as an expansion for {2MwY ^ which 
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we reproduce. One has 
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where A = 2^/^G';\a/7r), Tt = M'jj/{2Mtf, and Ty/ is defined in Section [TJ Here, a is 
Sommerfeld's fine-structure constant, Gp is Fermi's constant, A'c = 3 is the number of 
quark colours, and Qt = 2/3 is the electric charge of the top quark in units of the positron 
charge. The Higgs-boson mass entering the expansion parameters Tt and Ty/ partly arises 
from the couplings involving one Higgs and two Goldstone bosons, but also from the 
expansion in the external momenta due to the kinematic relation (gi + ^2)^ = ^'h- 

We wish to note that the result for the process H ^ gg is simply obtained by setting 
Ay) = and performing the substitution aNcQ^ — > asV^, where ag is the strong-coupling 
constant. 

In the remainder of this section, let us discuss the counterterm contributions needed 
for our two-loop analysis. We adopt the on-mass-shell scheme and regularize the ultra- 
violet divergences by means of dimensional regularization, with D = 4 — 2e space-time 
dimensions and 't Hooft mass scale /i. We use the anti-commuting definition of 75. 

As will become clearer in the next section, it is important to treat the tadpole contri- 
butions properly in our calculation. For this reason, in the following, we list the corre- 
sponding contributions separately and mark them by the superscript "tad" . Inspection of 
the one-loop diagrams reveals that, to 0{GfM^), we have to renormalize the Higgs-boson 
wave function and mass, the W^-boson mass, and the top-quark mass. The corresponding 
counterterms are defined through 
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Note that 6Zh is obtained from the derivative of the Higgs-boson self-energy and thus 
has no tadpole contribution. The mass counterterms are obtained from the corresponding 
two-point functions, where only the Mj-enhanced contributions have to be considered. Up 
to and including quadratic terms in Mf, we have 
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with xt = G F / {Stt'^ \/2) and A = 1/e — + ln(47r), where 7^; is Euler's constant. 
The H^-boson mass renormahzation is needed for the W, (p, and u propagators, where 
Mw enters as a parameter. Furthermore, also the HW^W^, H(j)'^W^, (p'^W^'j, and 
Hcp'^W^'y vertices contain Mw- The only vertex involving Mh is Hcp^cj)^, which induces 
two- loop contributions via 6Mh- Finally, Mt occurs in the top-quark propagator and in 
the Hti vertex. 

The wave function renormahzation and the renormahzation of a factor 1/Mw is com- 
mon to all one-loop diagrams. This allows for the definition of a universal factor ^3], 
which is finite in our calculation. It is given by 



The contributions of 0{GfM^) are obtained by considering all two- loop electroweak di- 
agrams involving a virtual top quark. This includes also the tadpole diagrams with a 
closed top-quark loop, which are proportional to M^. For arbitrary gauge parameter, 
this leads us to consider a total of 0(1000) diagrams. Some of them are depicted in 
Fig. 121 These diagrams naturally split into two classes. The first class consists of those 
diagrams where a neutral boson, i.e. a Higgs boson or a neutral Goldstone boson [x), is 
added to the one-loop top-quark diagrams. The exchange of a Z boson does not produce 
quadratic contributions in Mt. The application of the asymptotic-expansion technique to 
these diagrams leads to a simple Taylor expansion in the external momenta. 

This is different for the second class of diagrams, which, next to the top quark, also 
contain a W or (p boson and, as a consequence, also the bottom quark, which we take 




-NcXt. 
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3 Two-loop results 
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to be massless throughout the calculation. Due to the presence of cuts through light- 
particle lines, the asymptotic-expansion technique applied to these diagrams also yields 
nontrivial terms, as is exemplified in Fig. |H1 The first contribution on the right-hand side 
of the equation in Fig. IHl symbolizes the naive Taylor expansion in the external momenta. 
In the second contribution, the subdiagram to the right of the star has to be expanded 
in its external momenta, which also includes the loop momentum of the co-subgraph to 
the left of the star. The result of the expansion is inserted as an effective vertex, and 
the remaining integration over the second loop-momentum is performed after a further 
expansion in qi and q2- The latter is allowed, since we work in the limit {2Mw)'^ ^ Mfj. 
We wish to mention that such contributions develop terms, which cancel in the final 
result only in combination with the terms from the counterterms of Eq. (|7j) and the 
genuine two-loop tadpole diagrams. For this reason, it is crucial to include the latter 
in our calculation. As a further comment, we note that, unlike the example of Fig. 01 
it can also happen that the co-subgraph only involves massless bottom quarks, and no 
expansion in the external momenta is allowed. In our calculation, the contributions from 
such co-subgraphs vanish. 

As already mentioned above, we use in our calculation a general gauge parameter 
related to the W boson and verify that our final result is independent of C,w- We do this in 
the limit of large and small values of To this end, we apply the asymptotic-expansion 
technique in the following four limiting cases 

(t) > = ^wM^ > M|, 
(tt) > > ^wM^ > Ml, 
(m) > ^wM^ > > Ml, 

(iv) iwMl > M^ > M^ > M|, (9) 

where the inequality M^ ^ M| has to be understood in a formal sense, as, in practice, 
one has {2MwY ^ ^h- The result we obtain by asymptotic expansion is an expansion 
of the exact result for the 0{GfM^) contribution in powers of tw- In case (z), where 

= I5 which corresponds to 't Hooft-Feynman gauge, we are able to evaluate the first 
five terms of this expansion. In cases (u)-(if ), we compute the first two expansion terms 
and find agreement with the result obtained for C,w = 1- 

Our final result for A^ll emerges as the sum 

^« = ^« + A^Sx + -^wl, (10) 

where A^^ is the universal contribution induced by the one- loop term 6u of Eq. A^^]x 

is the two-loop contribution involving virtual H and x bosons, and A^)^ the remaining 

two- loop contribution involving virtual W and (p bosons. In A^^^^ and -Aj^^^, also the 
corresponding counterterm and tadpole contributions are included. For the individual 
pieces, we find 

(1) _ . /_329 _ 77 _ 19 2 _ 58 3 _ 1312 4 \ 
A^ -^A,Xi(^ ^^Tw ^^T^ 225^^ 7425^^ + "7' 
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Figure 2: Sample Feynman diagrams contributing at the two-loop electroweak order to 
the process H — > 77. 



.(1) 22 76 2 232 , 5248 4 \ 

where A is defined below Eq. dH} and the ellipses indicate terms of 0{t^). Notice that 
the leading OlGpMf) term of A^^\ '^ot accompanied by an expansion in tw, since the 
contributing diagrams do not involve virtual W or (p bosons. On the other hand, detailed 
inspection reveals that there is also no expansion in the parameter Mfj/{2Mz)'^, contrary 
to what might be expected at first sight. Inserting Eq. (fTT|) into Eq. (jTUl), we obtain our 
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Figure 3: Diagrammatic asymptotic expansion of a Feynman diagram that produces 
terms. 

final result 

.(1) A.r ^367 11 19 2 58 3 1312 ^ \ 
Al{l = ANcXt \ Tw H H H rl H . (12) 

The correction of 0{GfM^) to r(/7 — » 77) was also considered in Ref. ^Hl- The 
expression found in that reference disagrees with our result. One reason is probably that 
the authors of Ref. |1Q only considered virtual (f) bosons, but disregarded virtual W- 
bosons. However, our calculation shows that both types of charged bosons have to be 
taken into account at the same time in order to arrive at an ultraviolet-finite and gauge- 
parameter-independent result. In Ref. |i7], the dominant two-loop electroweak corrections 
to the Higgs-boson couplings to pairs of gauge bosons and hght fermions induced by a 
sequential isodoublet of ultraheavy quarks (A, B) were investigated by means of a low- 
energy theorem ITH] . In that paper, also a result for the 0{GfM^) correction to 
r{H 77) is presented, which is obtained by sending the mass of the fourth-generation 
down quark (Mb) to zero at the end of the calculation, which is performed assuming 
that Mb ^ Mw- This result also deviates from the one obtained above. Detailed 
inspection reveals that this difference may be attributed to the interchange of limits 
performed in Ref. ^7j, which is not justified in the case under consideration, although such 
a procedure is known to lead to correct results in simpler examples. In fact, reanalyzing 
the limiting case Ma ^ Mb ^ M]v using the asymptotic-expansion technique, we are 
able to reproduce the terms of Eqs. (48)-(54) in Ref. ^7] that survive in this limit and, 
at the same time, to identify contributions that do not occur if Mb = is imposed at the 
outset of the calculation. 

For completeness, we also specify our corresponding result for T{H ^ gg). Its evalu- 
ation is significantly simpler, since, at one loop, only the top-quark diagrams contribute. 
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Consequently, to obtain the leading two-loop correction proportional to GpM^, we only 
have to consider the exchange of H, x, and bosons. This only requires a naive Taylor 
expansion in qi and q2. We obtain 

-^s = ^9^* Q ' (13) 

with Agg = 23/4G]/^(a^/7r). This result is in agreement with Ref . |^ , where a low-energy 
theorem [TU HH] was used. 



4 Numerical results 

We are now in a position to present our numerical results and to assess the convergence 
properties of our expansions in and tw- We use the following numerical values for 
our input parameters ^: Gp = 1.16639 x 10"^ GeV-^, Mw = 80.423 GeV, and Mt = 
174.3 GeV. 

We first consider the one-loop amplitudes Af'^ and A^) induced by virtual top quarks 
and W bosons, respectively, for which exact results are available. They are shown in 
Figs, m and El as functions of Tt and tw, respectively. The solid curves indicate the 
exact results, while the dashed curves represent the sequences of approximations that 
are obtained by successively including higher powers of and t\y, respectively, in the 
expansions. The vertical lines encompass the intermediate-mass range of the Higgs boson, 
Mw < Mh < 2Mvy. In Fig. El the second vertical line coincides with the right edge of 
the frame. From Fig. EJ we observe that the approximation consisting of the first three 
terms of the expansion in Tt is practically indistinguishable from the exact result up to 
Tt ~ 0.25. The same is true for the sum of the first five expansion terms up to Tt ~ 0.5, 
which corresponds to Mh ~ 245 GeV. In the case convergence is slightly worse, 

since Mh = 2Mw corresponds to T]v = 1 and the exact result behaves like ^yl — tw in 
this limit. Nevertheless, for Mh = 120 GeV, 140 GeV, and 2Mwi the approximation by 
five expansion terms deviates from the exact result by as little as 0.3%, 1.6%, and 19.9%, 
respectively. 

Having demonstrated the fast convergence of the expansions in Tt and t\y of the one- 
loop amplitudes A^f^ and A^ , respectively, we now proceed to the two-loop electroweak 
amplitude ^j^y, for which we found the leading 0{GfM^) term together with its sub- 
leading mass corrections through (9(r^). This corresponds to a first approximation and 
four improvements, which are visualized by the five dashed curves in Fig. IHl As in Fig. El 
the dotted vertical line and the right edge of the frame enclose the intermediate-mass 
range of the Higgs boson. We again observe rapid convergence. The goodness of our best 
approximation may be assessed by considering its relative deviation from the second best 
one. For Mh = 120 GeV, 140 GeV, and 2Mw, this amounts to 0.3%, 1.0%, and 2.8%, 
respectively. The situation is very similar to the one in Fig. El In fact, the correspond- 
ing figures for A^) are 0.4%, 1.1%, and 3.1%. We thus expect that the goodness of the 
approximation of A^ll by the expansion through C(r^) is comparable to the case of A^) ■ 
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Figure 4: Top- quark- induced one-loop amplitude Af' normalized to AN^Ql as a function 
of Tt. The solid curve indicates the exact result, while the dashed curves represent the 
sequence of approximations that are obtained by successively including higher powers of 
Tt in the expansion. The dotted vertical lines encompass the intermediate-mass range of 
the Higgs boson. 

For the comparison with future measurements of T{H 77), all known corrections 
have to be included in Eq. Q. In this connection, it is interesting to compare the size 
of the new 0{GfM^) electroweak correction with the well-known 0{as) QCD one jl]. 
This is done in Fig. ((7j), where the respective corrections to T{H ^ 77) are displayed as 
functions of M^- As in Figs. ElandlHl the dotted vertical line and the right edge of the 
frame margin the intermediate-mass range of the Higgs boson. We observe that, within 
the latter, the O^GpM^) correction slightly exceeds the 0{as) one in magnitude, a rather 
surprising finding. Due to the sign difference, the two corrections practically compensate 
each other. The two-loop electroweak correction induced by light-fermion loops, which 
has become available recently , is also negative, but has a slightly smaller size than the 
0{GfM^) correction. 
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Figure 5: PF-boson-induced one-loop amplitude normalized to ^ as a function of 
Tw- The solid curve indicates the exact result, while the dashed curves represent the 
sequence of approximations that arc obtained by successively including higher powers of 
Tw in the expansion. The dotted vertical line and the right edge of the frame encompass 
the intermediate-mass range of the Higgs boson. 

5 Conclusions 

We calculated the dominant two-loop electroweak correction, of 0{GfM^), to the partial 
width of the decay into two photons of the SM Higgs boson in the intermediate mass 
range, < Mh < 2Mvi/, where this process is of great phenomenological relevance for 
the searches at hadron colliders of this elusive missing hnk of the SM. 

We evaluated the relevant Feynman diagrams by application of the asymptotic-expan- 
sion technique exploiting the mass hierarchy 2Mt ^ 2Mw ^ Mh- In this way, we 
obtained an expansion of the full 0{GfM^) result in the mass ratio tw — M'^/{2MwY 
through 0{t^). 

The convergence property of this expansion and our experience with the analogue 
expansion at the Born level, where the exact result is available for reference, lead us 
to believe that these five terms should provide a very good approximation to the exact 
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Figure 6: Top- quark- induced two-loop electroweak amplitude A[]1 normalized to ANcXt 
as a function of r^. The dashed curves represent the sequence of approximations that 
are obtained by successively including higher powers of in the expansion. The dotted 
vertical line and the right edge of the frame encompass the intermediate-mass range of 
the Higgs boson. 

result for Mh ^ 140 GeV. By the same token, the deviation of this approximation for the 
0{GfM^) amplitude A^^y from the unknown exact result for this quantity is likely to 
range from 2% to 20% as the value of Mh runs from 140 GeV to 2Mw- 

In the intermediate Higgs-boson mass range, the 0{GfM^) electroweak correction 
reduces the size of T{H — > 77) by approximately 2.5% and thus fully cancels the positive 
shift due to the well-known 0{as) QCD correction 

As a by-product of our analysis, we also recovered the 0{GpMf) correction to the 
partial width of the decay into two gluon jets of the intermediate-mass Higgs boson, in 
agreement with the literature |19j . 
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OiGpM^) electroweak correction (solid curve) is compared with the 0{as) QCD one 
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